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Abstract 



QCD in two dimensions is investigated using the improved fermionic lattice 
Hamiltonian proposed by Luo, Chen, Xu, and Jiang. We show that the 
improved theory leads to a significant reduction of the finite lattice spacing 
errors. The quark condensate and the mass of lightest quark and anti-quark 
bound state in the strong coupling phase (different from t'Hooft phase) are 
computed. We find agreement between our results and the analytical ones in 
the continuum. 
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I. INTRODUCTION 



Wilson's lattice formulation of QCD is the most reliable and powerful non-perturbative 
approach to the physics of strongly interacting particles, but its progress was hampered by 
systematical errors mainly due to the finite value of the lattice spacing a. 

Wilson's gluonic action differs from the continuum Yang-Mills action by order of 0(a 2 ), 
while the error of Wilson's quark action is bigger, i.e., the error being of the order 0(a). 
Due to quantum effects, there is an additional problem called "tadpole" . If a and the bare 
coupling constant g were small enough, these errors would be negligible. Unfortunately, even 
in recent 4-dimensional lattice QCD simulations on the most powerful computers, the lattice 
spacing a is larger than 0.1/m and the lattice coupling gi at is bigger than 0.9. For these 
lattice parameters, violation of scaling is still obvious and the extrapolation of the results 
to the a — > limit induces unknown systematic uncertainties when extracting continuum 
physics. 

One of the most efficient ways for reducing these systematic errors is the Symanzik im- 
provement program JlJ: adding local, nearest neighbor or next-nearest-neighbor interaction 
terms to the Wilson lattice theory so that the finite a errors become higher order in a. 
During recent years, application of the Symanzik program has been a major topic. There 
have been several proposals on this subject: 

(a) For the gauge sector, Lepage proposed a tadpole improved action 0, reducing the 
errors from 0(a 2 ) to 0(a 4 ). Luo, Guo, Kroger, and Schiitte constructed a tadpole improved 
Hamiltonian ||[|] with the same accuracy. 

(b) For the fermionic sector, Hamber and Wu (one author of the present paper) proposed the 
first improved action ||, by adding next-nearest-neighbor interaction terms to the Wilson 
quark action so that the 0(a) error is canceled. There have been some numerical simulations 
|6|-[8|] of hadron spectroscopy using the Hamber- Wu action. Sheikholeslami and Wohlert || , 
Liischer and Weisz [|l(J worked on fermionic improvement and recent progress has been 
achieved on non-perturbative improvement via implementation of PCAC by the ALPHA 



collaboration |TT| . Luo, Chen, Xu, and Jiang (two authors of the present paper) constructed 



an improved Hamiltonian which was tested in the Schwinger model (QED 2 ). There 



are other possibilities for improving Wilson's quark theory [13]. The purpose of this work 



is to demonstrate numerically that fermionic improvement works also in the Hamiltonian 
formulation for the case of QCD in two dimensions. 

Like Quantum Electrodynamics in 2 dimensions (QED2) also Quantum Chromodynamics 
in 2 dimensions (QCD 2 ) has the properties of chiral symmetry breaking, an anomaly and 
confinement. However, the latter has gluonic self interactions, which makes it more similar 
to QCD in 4 dimensions (QCD 4 ) than to the Schwinger model. In addition to having 
confinement and chiral-symmetry breaking, QCD2 has much richer spectrum because of the 
non-abelian gauge interactions between "gluons" and "quarks" . Early in 1974, 't Hooft |T4j 
did pioneering work on QCD 2 using the 1/N C expansion (where QCD corresponds to the 
gauge symmetry group SU(N C )). He showed in the limit N c — > 00 that planar diagrams 
are dominant and they can be summed up in a bootstrap equation. From that he obtained 
a meson-like mass spectrum. This model has been extensively studied in the large Nc limit 
and also in the chiral limit m q —>■ (where m q is the free quark mass). Zhitnitsky |16[ has 
pointed out that there two distinct phases: 
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(i) Nq — > oo firstly, and then m q — > afterwards, 

(ii) m q — > 0, and then N c — * oo. 

The first case corrersponds to g << m q , which is the weak coupling regime. It descibes the 
phase considered by t'Hooft, where the following relation holds, 

Nq — > oo, g 2 Nc = const. m q >> g ~ . (1.1) 

\fN c 

t'Hooft found a spectrum of states (numbered by n) given by 

n T 2 2 9 2 Nc n x 

7T 

Zhitnitsky [l6j has computed the gluon condensate and the quark condensate. The quark 



condensate is given by 



<W = -W^p (1-3) 



The second case corresponds to g » m q which is the strong coupling regime. In the strong 
coupling regime the observed spectrum is different |T7HiTfl. Steinhardt JL9| has computed 



the masses of the following particles: soliton (baryon), anti-soliton (anti-baryon) and soliton- 



antisoliton (baryon-antibaryon) bound state. Bhattacharya [18] has shown that in the chiral 
limit there are free fields of mass 



M- a ^. (M, 
Note that for Nc = 1 this coincides with the Boson mass in the Schwinger model. Grandou 



et al. [2D| have obtained quark condensate in the chiral limit (for arbitrary g) and obtained 

<^> - -gN 3 c /2 (1.5) 
which confirms Zhitnitsky's result obtained in the weak coupling regime. 

Of course, QCD2 is much simpler than QCD4. It has been used to mimic the properties 
of QCD4 such as vacuum structure, hadron scattering and decays, and charmonium picture. 
Unlike the massless Schwinger model, unfortunately, QCD2 is no longer exactly solvable. 

The first lattice study of QCD 2 was done by Hamer [^TJ using Kogut-Susskind fermions. 



He computed for SU(2) the mass spectrum in the t'Hooft phase. In 1991, Luo et al. p2 
performed another lattice field theory study of this model (SU(2) and SU(3)) using Wilson 
fermions. As is shown later, the results for lattice QCD 2 with Wilson quarks were found 
to be strongly dependent on the unphysical Wilson parameter r, the coefficient of the 0(a) 
error term. The purpose of this paper is to show that in the case of QCD 2 , the improved 
theory [|12| can significantly reduce these errors. 

The remaining part of the paper is organized as follows. In Sect. |H|, we review some 
features of the Hamiltonian approach as well as the improved Hamiltonian for quarks pro- 
posed in Ref. [12]. In Sect. |T|, the wave functions of the vacuum and the vector meson are 



constructed, and the relation between the continuum chiral condensate and lattice quark 
condensate is developed. The results for the quark condensate and the mass spectrum are 
presented in Sect. and discussions are presented in Sect. 0. 
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II. IMPROVED HAMILTONIAN FOR QUARKS 



Although numerical simulation in the Lagrangian formulation has become the main 
stream and a lot of progress has been made over the last two decades, there are areas 
where progress has been quite slow and new techniques should be explored: for example, 
computation of the S-matrix and cross sections, wave functions of vacuum, hadrons and 
glueballs, QCD at finite baryon density, or the computation of QCD structure functions 
in the region of small xb and Q 2 . In our opinion the Hamiltonian approach is a viable 
alternative p3|p4| and some very interesting results p5|- pl|1 have recently been obtained. 



Many workers in the Lagrangian formulation nowadays have followed ideas similar to the 
Hamiltonian approach (where the time is continuos, i.e., a t = 0) by considering anisotropic 
lattices with lattice spacings at << a s . The purpose is to improve the signal to noise ratio 
in the spectrum computation 0. 

In the last ten years, we have done a lot of work [p!^ , |32} -|37H on Hamiltonian lattice field 



theory, where the conventional Hamiltonian in the Wilson fermion sector is used: 



Hf = H m + H k + H r , 



H m = m^t/)(x)^( 



x 



H k = ^-Y. t f J ( x )^U{x,k)4j(x + k), 

la x,k 



H r = Y a Y3i x )^ x ) ~ $(x)U(x, k)i/j(x + k)), (2.1) 

where a is now the spacial lattice spacing a s , U(x, k) is the gauge link variable at site x in 
the direction k = ±j (j is the unit vector), and 7_j = —jj, H m , Hk, H r are respectively 
the mass term, kinetic term and Wilson term. The Wilson term (r ^ 0), proportional to 
O(ret), is introduced to avoid the fermion species doubling, with the price of explicit chiral- 
symmetry breaking even in the vanishing bare fermion mass limit. As discussed in Sec. |, 
the 0(a) error in Hf indeed leads to lattice artifacts if a or gi at is not small enough. 

Similar to the Hamber-Wu action ||, where some next-nearest-neighbor interaction terms 
are added to the Wilson action to cancel the 0(ra) error, we proposed a 0(a 2 ) improved 



Hamiltonian in Ref. [|12jl : 



jjimproved jj _j_ jjimproved _|_ jjimproved 



Hr r ° ved = ^Tp(x)lkU(x, k)1,{x + k) 

x,k 



+ Y a E ${x)lkU{x, 2k)*P(x + 2k), 

x,k 
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H improved = JL^ ^(x)lf>(x) 
x,k 



f _ 

-Ci— ip{x)U{x, k)ip(x + k) 



-c 2 — ^ijj{x)U{x,2k)ijj{x + 2k). (2.2) 

2 fl x,k 

Here U(x, 2k) = U(x, k)U(x + k, k) and the coefficients b\, fe 2 , c\ and c 2 are given by 

4 14 1 

h = - , b 2 = --, ci = -, c 2 = --. (2.3) 
3 o 3 3 

These coefficients are the same for any d+1 dimensions and gauge group. The results shown 
in this paper correspond to this set of parameters. However, the following set of paramters 

6 1 = l,6 a = 0,c 1 = -,c 2 = -^ (2.4) 

where only the Wilson term is improved, gives very similar results. With the absence of the 
0(ra) error, we expect that we can extract the continuum physics in a more reliable way. 

Lattice QCD 2 in the Hamiltonian approach has some nice features which simplify the 
computations considerably: 

(a) The magnetic interactions are absent and the tadpole factor Uq = 1 due to the fact that 
Up — Z7J = 1. Therefore, there is only color-electric energy term in the gluonic Hamiltonian: 
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2a 



with j — I and a — 1, N% — 1. 

(b) The quantum (tadpole) effects, if any, are highly suppressed in the fermionic sector 
as 0(gf at a) — > 0(a 3 ) and in the gluonic sector as 0(gf at a 2 ) — > 0(a 4 ). The reason is the 
super-renormalizability of 1+1 dimensional theories, where 

giat = ga. (2.6) 
All this means that in 1+1 dimensions, classical improvement of the fermionic Hamil- 



tonian is sufficient. In Ref. |12| , the mass spectrum of the Schwinger model was used to 
test the improved program. In the following sections, we will provide evidence in QCD 2 to 
support the efficiency of the improved Hamiltonian. The results for the quark condensate 
and mass spectrum further confirm our expectation. 
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III. VACUUM AND VECTOR PARTICLE STATES 



The vacuum wave function is constructed in the same way as in JT2 



where 



\Q) = exp(iS)\0), 



N s 

trun 



S = 9 n S n 



n=l 



Si = 4 1 { x )lkU (x, k)i/j(x + k) 



S 2 =iY, i/> j (x)lkU(x, 2k)ij(x + 2k) 

x.k 



S3 = ^ { x )lkU{x, 3k)t/j(x + 3k) 

x.k 



(3.1) 



(3.2) 

with 9 n determined by minimizing the vacuum energy. |0) is the bare vacuum defined by 



£(a:)|0) = t7(x)|0) = ^(x)|0) = 0. 



(3.3) 



Here £ and rf are respectively the up and down components of the ip field. Such a form of 
the fermionic vacuum (|3.2|) has also been discussed extensively in the literature [22 ,[$2 -37]. 

One is interested in computing the wave function of the lowest lying energy state which 
is the flavor- singlet vector meson | V), similar to the case of the Schwinger model. The wave 
function is created by a superposition of some operators V n with the given quantum numbers 



V = iX)V( a 07iV'( a 0> 



Vi — iy j ip(x)-yiU (x, k)ip(x + k), 
V 2 = i J2 $(x)nU(x, 2k)ip(x + 2k), 
V 3 = i 2 $(x)nU(x, 3k)tp{x + 3k), 
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(3.4) 



acting on the vacuum state \Q), i.e., 

\V)= £f A n [V n -(Q\V n \n}]\n). (3.5) 



n=0 



The criterion for choosing the truncation orders Nf run and N^. un is the convergence of the 
results. An estimate for the vector mass My is the lowest eigenvalue min(Ey) of the 
following equations 



NY 

trun 



Yl ( H n 2 m - E vU n2ni )A ni = 0, 



ni=0 



det\H v -E v U v \ = 0, (3.6) 

where the coefficients A ni are determined by solving the equations, and the matrix elements 
H n 2 m and u n 2 m are defined by 

U n 2 nx = (Vn 2 \Vm) L - (3.7) 

Here the superscript L means only the matrix elements proportional to the lattice size L are 
retained (those proportional to higher order of L do not contribute). Detailed discussions 
can be found in p8| . We can estimate the continuum My from 



Mv _ aMv^ (3 g) 

9 9lat 

if the right hand side is independent of gi at . 

In lattice field theory with Wilson fermions, chiral symmetry is explicitly broken even in 
the bare vanishing mass limit. Therefore, the fermion condensate {ipi/j) f ree is non-vanishing 
and should be subtracted in a way described in Refs. p2t^5ip6|^9| : 



(iplp)sub = - (4>1p)free, (3.9) 

where for one flavor 

4 1 dE n , , 

I I \ ^ L tree 



/n/a r 
dpa{[—(l — cicospa — C2Cos2pa)Y 
-n/a & 
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+ [ Smpa (b 1 + b 2 cospa)} 2 } 



r 

/{[— (1 — Cicospa — C2Cos2pa)} 



2 



+ [^(6 1 + W)] 2 }'/ 2 . (3.10) 
a 

From (ipip)sub, we can get an estimate of the continuum quark condensate (^V) cont 



9 9lat 



(3.11) 



if the right hand side does not depend on the bare coupling gi at . 

It is well known that spontaneous chiral-symmetry breaking originates from the axial 
anomaly and there is no Goldstone pion in quantum field theory in 1+1 dimensions. There- 
fore, for Wilson fermions, one can not fine-tune the (r, m) parameter space as in QCD 4 to 
reach the chiral limit. However, the chiral limit can be approximated as the m — > limit as 
long as the lattice spacing error is sufficiently small. For Wilson fermions, as discussed in 
Ref. PD]| , this is not well justified for finite gi at and a because of the 0{ra) error, but for the 
improved theory, as shown in Refs. [0,0) this would be a reasonably good approximation 



since the chiral-symmetry breaking term is much smaller, i.e., 0(ra 



2\ 



IV. RESULTS FOR QUARK CONDENSATE AND VECTOR MASS 



To increase the accuracy of the techniques described in Sect. |T|, we include higher order 
contributions in Eq. (|3.2|) than those in Ref. so that better convergence of the results 
can be obtained. 

Figure 1 and Figure 3 show {tpip) su bl 'giat as a function of l/gf at in 2-dimensional lattice 
SU(2) and SU(3), respectively, gauge theories with Wilson fermions. Figure 5 and Figure 7 
show a My I 'gi a t as a function of l/gf at in 2-dimensional lattice SU(2) and SU(3), respectively, 
gauge theories with Wilson fermions. As one can see, the results for r = 1 deviate obviously 
from those for r = 0.1, which is attributed to the 0(ra) error of the Wilson term. 

The corresponding results from the Hamiltonian with improvement are presented in 
Figure 2, Figure 4, Figure 6, and Figure 8. One observes that the differences between the 
results for r = 1 and r = are significantly reduced. Most impressively, the data for the 
quark condensate coincide each other. A similar r test has also been used in |15] for checking 
the efficiency of the improvement program. 

To get an idea about how QCD 2 behaves for large N c , we have computed the quark 
condensate (tpip) i n the chiral limit for Nq = 3,4,5,6. This is shown in Figure 9. We 
have compared our numerical results with the theoretical result by Zhitnitsky ||16| , Eq. (|1.3j ). 
Although Zhitnitsky's result was obtained in the weak coupling phase, it qualitatively agrees 
with the strong coupling result of Ref. [|20|| . Remarkably, our results agree very well with 
Zhitnitsky's weak coupling result. Secondly, we show in Figure 10 the lowest lying mass 
of the mass spectrum again for Nc = 3,4,5,6 in the chiral limit. Note that this particle 
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corresponds to the vector particle in the Schwinger model. In QCD2 it corresponds to 
free particles (the mass of which remains finite when m q goes to zero) distinct from "sine- 
Gordon" soliton particles (the mass of which goes to zero when m q goes to zero, see |fL9|j). 
We compare our numerical result with the analytical strong coupling result by Bhattacharya 



T8|, Eq.(|1.4|). Again we find agreement. 



V. CONCLUSIONS 

In this work we have shown that fermionic improvement works also in the Hamiltonian 
lattice formulation for the case of QCD\ + i. We have computed the quark condensate and 
the flavor singlet vector mass using the 0(a 2 ) improved Hamiltonian for quarks [12|] proposed 



by Luo, Chen, Xu, and Jiang. In comparison with the results corresponding to the Wilson 
fermions without improvement, we indeed observe significant reduction of the finite lattice 
spacing error 0(ra). By comparison with analytical results for the quark condensate and 
the vector mass we find good agreement. In our opinion, this is the first lattice study of 
QCD\ + \ which gives results in the strong coupling phase (in contrast to the t'Hooft phase). 
In particular we present results for different gauge groups (Nq = 2,3,4,5,6). We believe 
that the application of the Symanzik improvement program to QCD in 3+1 dimensions will 
be very promising. 
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FIGURES 
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FIG. 1. x = —(ip^svb/idlatNc) versus l/gf at for Nc = 2 with Wilson fermions. Crosses: 
r = 0.1, Diamonds: r = 1. 
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FIG. 2. x = —(ipi())sub/(giatN c ) versus l/sf at for Nc = 2 with improved Wilson fermions. 
Crosses: r = 0.1, Diamonds: r = 1. 
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FIG. 3. x = —(^^sub/idlatNc) versus l/.gf at for iVc = 3 with Wilson fermions. Crosses: 
r = 0.1, Diamonds: r = 1. 




FIG. 4. % = —(ip^}sub/(9iatN c ) versus l/sf at for iV"c = 3 with improved Wilson fermions. 
Crosses: r = 0.1, Diamonds: r = 1. 
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FIG. 5. aMy / giat versus l/gf at for Nc = 2 with Wilson fermions. Crosses: r = 0.1, Diamonds: 
r = 1. 
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FIG. 6. aMy/giat versus l/gf at for Nc = 2 with improved Wilson fermions. Crosses: r = 0.1, 
Diamonds: r = 1. 
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FIG. 7. aMy/giat versus l/gf at for iVc = 3 with Wilson fermions. Crosses: r = 0.1, Diamonds: 
r = 1. 
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FIG. 8. aMy/gi a t versus l/sf at for iVc = 3 with improved Wilson fermions. Crosses: r = 0.1, 
Diamonds: r = 1. 
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FIG. 9. p = —(tp^cont/igN^/ 2 ) in the continuum versus The error bars are estimated 



from the data for different r. The full line gives the analytical result from Ref. [16|. 




FIG. 10. n = M v /(gyJN c + 1) in the continuum versus 1/Nq. The error bars are estimated 
from the data for different r. The full line gives the analytical result from Ref. [18|. 
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